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1.  Contract  Description. 

V 

This  project  is  designed  to  investigate  a  wide  range  of 
phenomena  related  to  energetic  materials  in  their  solid  phase  and 
also  on  other  molecular  solids  which  serve  as  prototypes.  These 
studies  are  entirely  theoretical  in  nature  and  are  employing 
exclusively  nan-parametric  ab-initio  techniques. 

In  the  case  of  solid  fuels  or  explosives  the  area  of  prime 
concern  is  the  initiation  of  combustion  or  detonation.  This  may 
be  accomplished  by  means  of  shock,  or  of  high  temperature  or 
both. Thus  we  need  to  be  able  to  study  these  systems  under  extreme 
conditions  of  temperature  or  pressure.  In  order  to  accomplish 
this  objective  it  is  necessary  to  extend  the  normal  solid  state 
methods  employed  for  studying  the  electronic  or  vibronie 
structures  as  such  theories  are  essentially  T«0  theories. 

If  one  is  able  to  explain  the  initiation  process,  then  it  is 
essential  to  be  able  to  explain  how  one  might  sustain  the 
compustion  as  well.  One  likely  mode  of  sustaining  the  reaction  is 
an  auto-catalytic  mechanism. 

Finally  one  need  examine  those  properties  of  the  solid  which 
contribute  to  the  stability  of  the  system,  and  if  possible  how 
such  properties  ara  inf luenced/control led  by  impurities  or 
defects  or  their  aggregates. 


2.  Partial  List  of  Unsolved  Questions  of  Scientific  Interest. 

If  one  is  to  except  the  solid  rare  gases  and  such  inherently 
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one-dimensional  systems  as  polymers,  very  few  theoretical  studies 
have  been  performed  on  crystaline  molecular  solids.  The  majority 
of  the  extant  studies  are  by  means  of  cluster  simulations  rather 
than  by  use  of  extended  system  techniques (such  as  energy  band 
theory).  In  addition  the  majority  of  studies  are  of  a  partly  (or 
even  totally)  empirical  nature.  The  one  current  exception  to  this 
is  the  case  of  solid  CH4,  which  has  been  studied  in  a  Hartree— 
Fock  model  some  10  years  ago.  Unfortunately  the  calculation  made 
some  unacceptable  approximations  to  the  Fock  matrix  elements  and 
is  not  wholly  reliable.  This  and  other  criticisms  of  the  previous 
work  are  detailed  in  the  accompaning  manuscripton  solid  CH4.  In 
addition  the  present  results  extend  past  knowledge  by  including 
explicit  treatment  of  correlation  effects  and  also  in  that  it 
includes  the  effects  of  hydrostatic  pressure. 

The  fundamental  theoretical  knowledge  of  even  such  a  simple 
energetic  solid  as  NO  is  to  all  intents  and  purposes  zero. 
Clearly  then  one  does  not  expect  (nor  does  one  find)  any  detailed 
theoretical  understanding  of  the  fundamental  processes  in  more 
complicated  and  more  important  (technologically)  solid  energetic 
systems. 


3.  Scientific  Approach  and  Progress  in  its  Implementation 
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Two  fundamental  and  complementary  approaches  are  to  be 
employed  in  this  study.  The  first  of  these  is  applicable  to 
systems  having  translational  invariance  is  energy  band  theory, 
the  second  is  more  applicable  to  those  systems  lacking  even  one0r 


di Mnii on  with  translational  invariancs  such  as  systems  with 
point  defects  or  impurities,  this  bsing  extended  clustsr 
nodsling.In  this  latter  approach  ths  extended  system  is  modeled 
by  a  group  of  molecules  , twenty  or  so  in  number,  which  are 
bounded  by  a  self -consistently  determined  potential  to  simulate 
the  remainder  of  the  system.  This  cluster  is  then  solved  self- 
consistently  in  the  Unrestricted-Hartree-Fock  (UHF)  approximation 
separately  for  each  state  in  question.  The  virtual  UHF  orbitals 
are  retained  and  employed  as  a  basis  set  for  computing  the 
correlation  energy  of  the  states  in  question  by  means  of  Many- 
Body-Perturbation-Theory  <MBPT> .  This  is  not  an  unusual  method 
for  studies  of  molecular  ground  states,  but  what  is  unusual  here 
is  its  employment  on  such  a  large  system  and  also  for  excited 
states.  In  order  to  develop  the  excited  state  technology  we  have 
been  looking  at  some  open  shell  singlets.  In  so  doing  it  becomes 
necessary  to  employ  several  projection  techniques  to  prevent 
variational  collapse.  This  is  scientifically  important  but  not 
yet  critical  to  this  program  sothat  this  particular  study  is 
given  very  low  priority.  None  the  less,  early  results  on  some 
atoms  and  small  molecules  indicate  an  ability  to  compute  the 
energy  of  excitation  from  ground  state  singlets  to  excited  open 
shell  singlets  to  an  accuracy  of  better  than  0.2  eV. 

Several  useful  studies  on  CH4  solid  have  been  performed 
using  these  techniques.  These  include  a  study  of  how  the  C-H  bond 
length  varies  as  a  function  of  the  C-C  bond  distance.  That  is, 
the  internal  rearrangements  of  the  CH4  molecule  in  the  solid  are 
studied  as  a  function  of  hydrostatic  pressure.  This  calculation 


permits  us  to  deduce  several  vary  useful  parameters.  Thasa 
includa  tha  effective  CH4-CH4  force  law,  the  C-H  force  law  and 
tha  H-H  forca  law.  To  a  tolarabla  approxiaation  all  thasa  forcas 
ara  harmonic  and  tha  assantial  CH4-CH4  forca  constant  is  quits  a 
bit  softar  than  any  of  tha  othars.  This  is  usaful  in  that  it 
parmits  one  to  partition  tha  vibrational  modes  into 
inter molecular  and  intramolecular  modes.  This  has  been  dona  and 
tha  vibronic  properties  of  the  bulk  system  deduced.  Defects  were 
introduced  in  the  bulk  and  on  the  surface  by  removing  a  CH4 
molecule,  and  the  vibrational  energies  recomputed.  These  results 
ara  interesting  and  we  believe  significant.  We  find  that  a  bulk 
vacancy  introduces  a  local  soft  mode  into  the  vibrational 
structure,  whereas  the  surface  vacancy  does  not.  Therefore  we 
conclude  that  the  bulk  vacancy  may  serve  as  an  energy  trap  but 
tha  surface  one  is  not  likely  to  do  so.  Clearly  this  point  needs 
further  study  on  several  more  interesting  cases,  but  this 
auggnts  that  the  surface  vacancy  is  not  likely  to  act  as  a  point 
for  initiation  of  combustion  if  tha  enargy  transfer  mode  is  the 
vibrational  spectrum.  Clearly  we  will  pursue  this  study  further 
and  extend  it  to  NO  as  soon  as  practical.  It  is  intended  to  write 
up  the  basic  results  of  tha  cluster  studies  for  CH4  for 
publication  in  the  near  future. 

The  implementation  of  ab  initio  energy  band  methods  to  tha 
study  of  molecular  solids  seems  to  at  last  be  a  usaful  success. 
Tha  method  of  choice  is  tha  linear-combination-of-molecular- 
orbitals  (LCNO)  approach  for  technical  reasons  detailed  in  the 
Pri w  reports.  The  initial  attempt  to  implement  this  idea  was 
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technically  a  success  ,but  -far  too  slow  in  practice  to  cope  with 
soiids  more  complex  than  A r,  or  CH4.  Therefore  with  the  goals  of 
this  project  in  mind,  I  determined  to  reformulate  the  entire  LCMO 
approach  to  gain  the  necessary  speed.  This  has  been  successfully 
completed,  and  a  code  now  exists  which  should  easily  allow  one  to 
study  energy  bands  in  systems  as  complex  as  NO  or  even  more 
complex.  The  energy  band  studies  here  are  at  the  HF  level.  If  one 
is  to  explain  experiment  or  have  quantitative  predictive  ability, 
the  inclusion  of  correlation  effects  is  necessary.  There  are  two 
essential  problems  here,  the  first  is  that  the  correlation  energy 
for  a  solid  is  infinite, but  has  finite  differences  upon 
excitation.  Thus  any  attempt  to  evaluate  the  correlation  energies 
in  two  states  is  foredoomed  to  failure,  and  one  must  develop 
techniques  to  evaluate  the  differences  directly.  The  second 
problem  is  that  of  size  consistency.  This  problem  is  somewhat 
avoided  as  is  seen  in  the  preliminary  manuscript  on  solid  CH4, 
when  one  computes  changes  in  energy  rather  than  the  energies 
directly.  The  changes  are  at  least  in  part  not  extensive 
quantities  as  we  show.  Those  parts  which  are  extensive  are  then 
evaluated  by  means  of  MBPT  methods.  The  essential  interesting 
result  for  solid  CH4  is  that  the  energy  band  gap  is  a  highly 
sensitive  function  of  the  lattice  constant.  The  lowest  gap  is 
indirect  and  has  a  value  of  13.17  eV  at  normal  lattice  constant 
and  decreases  precipitously  to  10.86  eV  when  the  lattice  constant 
is  decreased  by  5%.  This  suggests  for  narrow  gap  solids  that 
pressure  plus  elevated  temperature  could  excite  quite  a  few 
molecules  out  of  their  electronic  ground  state.  The  utility  of 
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this  is  further  sssn  if  ms  include  the  electron-hols  interaction 
betMeen  the  excited  electron  and  the  hole  left  in  the  valence 
band.  This  interaction  causes  the  two  entities  to  fore  an  exciton 
and  this  drives  a  local  distortion  of  the  CH4  molecule  from  the 
ideal  Td  geometry  to  D2h.  In  this  change  one  enables  the 
separation  of  the  CH4  into  two  reactive  elements,  H2  and  CH2,  if 
one  further  perturbs  the  system.  This  suggests  that  electronic 
excitations  may  play  a  significant  role  in  the  oxidation  of  such 
systems.  The  full  details  of  this  study  are  attached.  It  is 
expected  to  pretty  up  a  few  details  of  the  codes  and  then  start 
on  solid  NO. 
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